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Abstract 



We present results from a search for anomalous WW and WZ production 
in pp collisions at t/s = 1.8 TeV. We used pp — > eujjX events observed 
during the 1992-1993 run of the Fermilab Tevatron collider, corresponding to 
an integrated luminosity of 13. 7 ±0.7 pb _1 . A fit to the transverse momentum 
spectrum of the W boson yields direct limits on the CP-conserving anomalous 
WW^f and WWZ coupling parameters of —0.9 < An < 1.1 (with A = 0) and 
—0.6 < A < 0.7 (with Ak = 0) at the 95% confidence level, for a form factor 
scale A = 1.5 TeV, assuming that the WW^y and WWZ coupling parameters 
are equal. 

PACS numbers: 14.70.Fm, 13.40.Em, 13.40.Gp 
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Self-interactions of the electroweak gauge bosons such as WW^y and WW Z are a direct 
consequence of the non-Abelian gauge symmetry of the Standard Model (SM). The strength 
of these interactions can be directly measured by studying gauge boson pair production in pp 
collisions 0. For example, WW production is sensitive to the WWj and WW Z couplings, 
and WZ production is sensitive to the WWZ coupling. Any deviation of these couplings 
from their SM values indicates physics beyond the SM. 

The WW 7 and WWZ interactions are generally described by a Lagrangian with fourteen 
independent coupling parameters 0, of which k 7 and A 7 for the WWj vertex and k z and 
\z for the WWZ vertex are usually studied. These are CP-conserving coupling parameters. 
In the SM, Ak 7 (= k 7 — 1) = A 7 = Ak z (= k z — 1) = Xz = 0, and the production cross 
section for pp -> W + W~X (W ± ZX) at Vs = 1.8 TeV is 9.5 (2.5) pb §. 

Non-SM (i.e., anomalous) couplings dramatically increase the production cross section, 
and enhance the transverse momentum spectrum of the W boson (pr) for large values of 
p^. Therefore, a study of the p^ spectrum of WW (W Z) production leads to a sensitive 
test of the size of the WWj and WWZ couplings. With anomalous couplings, some helicity 
amplitudes of the pp — > WW(WZ) processes grow with s, the square of the invariant mass 
of the WW {W Z) system, and cause the cross section eventually to violate tree level S- 
matrix unitarity. To avoid this, the anomalous couplings are commonly parametrized as 
dipole form factors with a cut-off scale A: Ak(s) = A«/(l + s/A 2 ) 2 , A(s) = A/(l + s/A 2 ) 2 . 
Consequently, limits on the couplings An and A are dependent on the choice of A. 

The D0 collaboration has reported limits on WW'y and WWZ anomalous couplings 
from two processes using data from the 1992-1993 Fermilab Tevatron collider run with pp 
collisions at ^/s = 1.8 TeV: on the WW'y coupling from a measurement of W'y production || 
and on the WW'y and WWZ couplings from a search for W boson pair production in the 
dilepton decay modes ||. In this letter we present a new, independent, determination of 
limits on the WW'y and WWZ anomalous couplings obtained from a search for pp — ► WWX 
followed by W — > ev and W —>■ jj, and pp — > WZX followed by W — » ev and Z — > jj, 
where j represents a jet. Owing to the limited jet energy resolution, we cannot distinguish 
WZ events from WW events. This analysis uses the same data set as in the papers cited 
above, corresponding to an integrated luminosity of 13.7 ±0.7 pb -1 . The CDF collaboration 
has reported a similar measurement ||. 

The D0 detector and data collection systems are described in Ref . . The basic elements 
of the trigger and reconstruction algorithms for jets, electrons and neutrinos are given in 
Ref. [§. 

The WW, WZ — > eujj candidates were selected by searching for events containing a 
W — > ev decay and at least two jets consistent with W — > jj or Z — > jj. The data sample 
was obtained with a trigger which required an isolated electromagnetic (EM) calorimeter 
cluster with transverse energy Ex > 20 GeV. In offline event selection, this EM cluster 
was required to be within \r]\ < 1.1 in the central calorimeter, or 1.5 < \rj\ < 2.5 in the 
end calorimeters, where rj is the pseudorapidity, defined as rj = — ln[tan6 l /2] and 9 is the 
polar angle with respect to the proton beam direction. Such an EM cluster was identified 
as an electron if (i) the ratio of EM energy to the total shower energy was greater than 0.9; 

(ii) the lateral and longitudinal shower shapes were consistent with those of an electron; 

(iii) the isolation variable of the cluster was less than 0.1, where isolation is defined as 
/ = (E tot (0A) — E EM (0.2)) / E EM (0.2) , and E tot (0A) is the total calorimeter energy inside 
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a cone of radius 7Z = J (A?]) 2 + (A0) 2 = 0.4, where is the azimuthal angle around the 
beam axis and _Eem(0.2) is the EM energy inside a cone of radius 0.2; and (iv) a matching 
track was found in the drift chambers. The W — > ev decay was identified by an electron 
with Ej, > 25 GeV and missing transverse energy E? > 25 GeV forming a transverse mass 

Mf? = [2E^E T (1 - cos0 e!y )/c 4 ] 1/2 > 40 GeV/c 2 , where (p eu is the angle between the E e T and 
$t vectors. 

Jets were reconstructed using a cone algorithm with radius 1Z = 0.3, and were required 
to be within \rj\ < 2.5. The jet energies were corrected for detector effects: jet energy 
scale calibration and out-of-cone showering; for energy from the underlying event; and for 
energy loss due to out-of-cone gluon radiation |7j. We required that a candidate event 
contain at least two jets with E 3 T > 20 GeV and that the dijet invariant mass (the largest 
invariant mass if there were more than two jets with E° T > 20 GeV in the event) satisfy 
50 GeV/c 2 < rrijj < 110 GeV/c 2 , consistent with W and Z boson masses. Monte Carlo 
studies showed that the standard deviation of the dijet invariant mass distribution of the 
signal events is 15 GeV/c 2 . The above selection criteria yielded 84 candidate events. 

The trigger and electron selection efficiencies were measured || using Z — > ee events. 
The product of these efficiencies was found to be 0.78 ± 0.02 in the central calorimeter and 
0.62 ± 0.01 in the end calorimeters. The W — > jj selection efficiency was parametrized as 
a function of p^ , as shown in Fig. p]. This was estimated using events generated with the 
ISAJET |J and pythia ]TUj programs, followed by a detailed simulation of the D0 detector 
based on the geant package |1]J and application of our event selection criteria. The roll-off 
of efficiency in the p^ region beyond ~ 350 GeV/c is due to merging of the two jets from 
the W or Z boson. The use of a cone size as narrow as TZ = 0.3 for jet reconstruction 
ensures that the loss of efficiency occurs only for W and Z bosons with transverse momenta 
greater than expected from the coupling parameter values studied. The Z — > jj efficiency 
was obtained in a similar manner. In estimating the detection efficiencies of the WW(W Z) 
process, we used the W(Z) — > jj efficiency obtained from isajet, which is smaller than 
that from pythia and therefore gives more conservative limits on the coupling parameters. 

We calculated the overall event selection efficiency as a function of the coupling param- 
eters using a fast detector simulation program which incorporates the efficiencies described 
above and the detector resolutions [Jj]]. The WW(WZ) events were generated with the 
Monte Carlo program of Zeppenfeld [1,12], in which the processes were generated to leading 
order, and higher order QCD effects on the cross section were approximated by a factor 
K — 1 + §7ra: s = 1.34. We included the px distribution of the isajet WW events in the 
simulation of the WW(WZ) production. We calculated the total efficiency for SM couplings 
to be 0.15 ±0.02 for WW and 0.16 ±0.02 for WZ. The error is 13%, which is the addition in 
quadrature of the uncertainties on the electron trigger and selection efficiencies (2%), on the 
W(Z) — > jj efficiency due to the difference between the isajet and pythia programs (9%), 
the statistics of the Monte Carlo samples (4%), Et smearing (6%) and jet energy scale (6%). 
Therefore, 3.2 ± 0.6 WW and WZ events are expected based on the SM (2.8 ± 0.6 WW 
events plus 0.4 ±0.1 W Z events). The error here is the sum in quadrature of the uncertainty 
in the efficiency above and that of the higher order QCD corrections to the signal prediction 
(14%) @. 

The background estimate, summarized in Table I, includes contributions from: multijet 
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events, where a jet was misidentified as an electron and there was significant (mismeasured) 
missing transverse energy; W+ > 2j events with W — > eu; ti — > W + W~bb — > eujjX; 
WW (W Z) with W — > tv followed by r — > ei/9; and ZX — > eeX where one electron was 
not identified. 

The multijet background was estimated from the data by measuring the $t distribution 
of a background-dominated sample, which was obtained by selecting events containing an 
EM cluster that failed at least one of the electron identification requirements (ii) to (iv) 
described previously (shower shape, isolation or track- match). The $t distribution of this 
sample was scaled to match the candidate sample in the region GeV < $r < 15 GeV 
(before the $t requirement was applied), where the contribution of signal events is negligible; 
then the portion of this distribution passing the Ifir requirement {fir > 25 GeV) was taken 
as our estimate of the multijet background, giving 12.2 ± 2.3(stat.) ± l.l(syst.) events. The 



W+ > 2j background was estimated using the vecbos Monte Carlo program |13[, with 
Q 2 = m^, followed by parton fragmentation using the is A jet program and a detailed 
detector simulation. We normalized the number of VECBOS W+ > 2j events to the number 
of observed W+ > 2j events (after multijet events were subtracted) outside of the dijet 
mass signal region. This yielded 62.2 ±8.2(stat.) ± 10.1(syst.) W+ > 2j background events 
inside the rrijj signal region, where the statistical uncertainty is due to the size of the 
VECBOS W+ > 2j event sample and of all samples used to calculate the normalization 
factor (13%), and the systematic error is due to the normalization and to the jet energy scale 
correction (16%). The W+ > 2j cross section obtained with this normalization procedure 
was consistent with the VECBOS prediction. The backgrounds due to ti — > W + W~bb HMj , 
WW (W Z) — > rujj [l[j and ZX — > eeX were estimated using the isajet program followed 
by the geant detector simulation and found to be small. The total background from all 
sources was estimated to be 75.5 ± 13.3 events. Therefore, we observed no statistically 
significant signal above the background. Figure || shows the pr distribution of the ev system 
of the data events, background estimates, and Monte Carlo predictions of WW and WZ 
production for SM couplings and for one example with anomalous couplings. 

Using the efficiencies for SM WW and WZ production, the background-subtracted 
signal, the branching fractions B(W —>■ eu), B(W — > hadrons) and B[Z — > hadrons) 



from Ref. fL5 |, and assuming the SM ratio of the cross sections of pp — > W + W X and 



pp — > W ± ZX, we set an upper limit at the 95% confidence level (CL) on the cross section 
a{pp -> W+W~X) of 183 pb. 

The absence of an excess of events with high p^ excludes large deviations from the SM 
couplings. To set limits on the anomalous coupling parameters, we performed a binned 
likelihood fit on the entire p^ spectrum. For each pYp bin in the likelihood fit, and for a 
given set of anomalous coupling parameter values, we calculated the probability for the sum 
of the background and the predicted signal to fluctuate to the observed number of events. 
The uncertainties in the efficiency, background estimates, integrated luminosity and higher 
order QCD corrections to the signal prediction were convoluted in the likelihood function 
with Gaussian distributions. Our analysis is sensitive to anomalous couplings at both large 
and small s, since we do not require high p^, and therefore we retain events with small 



s [0. 

We obtained limits on the coupling parameters using four different assumptions of how 
the parameters are related to each other (a-d, below). The 95% CL contour limits are shown 
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as the inner curves in Fig. |3], along with the S-matrix unitarity limits, shown as the outer 
curves, which are obtained by evaluating all (i.e., W7, WW and WZ) processes. The 95% 
CL limits on the axes are listed in Table II. Two A values are used, A = 1.5 TeV for (a), (b) 
and (c), and A = 1.0 TeV for (d). These are values used in other measurements [3,5]. For 
(d), the results obtained with A = 1.5 TeV violate the S-matrix unitarity limit (not shown). 

In conclusion, we have searched for anomalous WW and WZ production in the eujjX 
decay mode, and set limits on the WWj and WWZ anomalous coupling parameters An 
and A. The limits on A are comparable to those measured using W'j production, whereas 
those on Ak from this analysis are significantly better. The results with assumption (c) 
are unique to WW {W Z) production since the WWZ couplings are not accessible with W^ 
production. The limits on both Ak, and A are significantly tighter than those from the 
analysis using WW to dilepton decay, due to the additional WZ production mode here, and 
to the larger branching fractions of W(Z) decaying to hadrons than that of W to ev or fiv 
in the dilepton analysis. 
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Science and Technology Policy (Russia), CNPq (Brazil), Departments of Atomic Energy 
and Science and Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry 
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TABLES 





TABLE I. Summary of evjj data 


and backgrounds. 






Background source 
multijets 

y v "t ^_ £j 

ti(m t = 180 GeV/c 2 ) 
VL1L,ITZ -► Ti/jj 
ZX -> eeA 
Total Background 


evjj events 
12.2 ±2.6 

DZ.Z ± lo.U 

0.87 ±0.12 
0.22 ±0.02 

o oo +0 - 34 

u - uu -o.oo 
75.5 ± 13.3 






Data 


84 






SM + 1LZ prediction 


3.2 ±0.6 




TABLE II. 95% CL limits with various assumptions. 


Assumptions 


A (TeV) 


Limits on axes 




(a) Ak 7 = Akz 
A 7 = Xz 


1.5 


-0.9 <Ak< 1.1 (A = 
-0.6 < A < 0.7 (A/e = 


= 0) 
= 0) 


(b) HISZ Q 


1.5 


-1.0 < Aft 7 < 1.3 (A 7 
-0.6 < A 7 < 0.7 (Ak 7 


= 0) 
= 0) 


(c) SM WWj 


1.5 


-1.1 < Ak z < 1.3 (Xz 
-0.7 < X z < 0.7 (Ak 2 


= 0) 
= 0) 


(d) SM WWZ 


1.0 


-2.8 < Aft 7 < 3.3 (A 7 
-2.5 < A 7 < 2.6 (Ak 7 


= 0) 
= 0) 



9 



FIGURES 



0.8 

£q 0.7 

T 0.6 

^ 0.5 

>^ 0.4 
o 

C 03 
'Q 0.2 









+ + ++ 






i 




^4 



50 100 150 200 250 300 350 

P T (W^jj) (GeV/c) 



400 450 500 



FIG. 1. Total efficiency for W — > jj selection as a function of p^, estimated using the isajet 
(solid) and pythia (dashed) generators followed by a full detector simulation. 




p T (W^ev) (GeV/c) 

FIG. 2. px distributions of the eu system: data (points), total background (solid), W+ > 2 jets 
background (dotted), and Monte Carlo predictions of WW and WZ production with SM (dashed) 
and non-SM (dot-dashed, Akz = Ak 7 = 2, Xz = A 7 = 1.5) couplings. 
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FIG. 3. Contour limits on anomalous coupling parameters at the 95% CL (inner curves) and 
limits from S-matrix unitarity (outer curves), assuming (a) Ak = Ak 7 = Akz, A = A 7 = Xz; 
(b) HISZ relations [17]; (c) SM WW-f couplings, and (d) SM WW Z couplings. A = 1.5 TeV is 
used for (a), (b) and (c); A = 1.0 TeV is used for (d). The SM prediction is Ak = 0, A = 0. 
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